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A frequentist asymptotic expansion method for error estimation is employed for a network of 
gravitational wave detectors to assess the amount of information that can be extracted from gravi¬ 
tational wave observations. Mathematically we derive lower bounds in the errors that any parameter 
estimator will have in the absence of prior knowledge to distinguish between the post-Einsteinian 
(ppE) description of coalescing binary systems and that of general relativity. When such errors 
are smaller than the parameter value, there is possibility to detect these violations from GR. A 
parameter space with inclusion of dominant dephasing ppE parameters (/3, b) is used for a study 
of first- and second-order (co)variance expansions, focusing on the inspiral stage of a nonspinning 
binary system of zero eccentricity detectible through Adv. LIGO and Adv. Virgo. Our procedure 
is an improvement of the Cramer-Rao Lower Bound. When Bayesian errors are lower than our 
bound it means that they depend critically on the priors. The analysis indicates the possibility 
of constraining deviations from GR in inspiral SNR (p ~ 15 — 17) regimes that are achievable in 
upcoming scientific runs (GW150914 had an inspiral SNR ~ 12). The errors on j3 also increase 
errors of other parameters such as the chirp mass AA and symmetric mass ratio p. Application is 
done to existing alternative theories of gravity, which include modified dispersion relation of the 
waveform, non-spinning models of quadratic modified gravity, and dipole gravitational radiation 
(i.e., Brans-Dicke type) modifications. 

PACS numbers: 04.30.Db, 04.80.Cc 


I. INTRODUCTION 

The advanced generation of the LIGO-Virgo network 
of interferometers m started collecting data in Septem¬ 
ber 2015 and provided the first detection of gravitational 
waves (GWs) [4], allowing to start testing General Rela¬ 
tivity (GR) beyond current constraints into strongly 
relativistic regimes ISHHl- In this paper we quantify the 
capability of laser interferometers to detect violations of 
GR, with a single detection of a compact binary coales¬ 
cence signal, by assessing if the minimal error on the pa¬ 
rameterized post-Einsteinian (ppE) parameters are larger 
than the separation of modified gravity values with re¬ 
spect to standard GR values. Error bounds are com¬ 
puted with the most accurate frequentist approach to 
date by computing the errors as inverse power series in 
the signal-to-noise ratio (SNR), where the first order is 
the inverse of the Fisher information matrix piTT]. In 
this paper we model GR violations with the ppE frame¬ 
work [aiiMis], which produces parametrized extensions 
of GR GW signals for the inspiral phase only of a binary 
compact coalescence in the absence of spin (similar ex¬ 
tensions are currently not available for the merger and 
ringdown phase as well as in the presence of spin). 

The square root of the inverse Fisher matrix diagonal 
elements, also known as the Cramer-Rao Lower Bound 
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(CRLB), is a lower limit in the error of any unbiased 
estimator in the absence of prior knowledge. In this re¬ 
gard the CRLB is a statement about the amount of in¬ 
formation available in the data regardless of the specihc 
parameter estimation scheme. There is however no guar¬ 
antee that any estimator is capable to actually attain the 
CRLB for part or the whole range of values the physical 
parameters can assume. Also, the CRLB only takes into 
account the curvature of the probability distribution of 
the data around the true value of the parameters and 
therefore does not include the role of secondary maxima 
in the calculation of the variance or mean square error of 
the estimators. The improved bound adopted here (based 
on second order asymptotics) is larger than the inverse 
Fisher matrices, known to underestimate errors in low- 
SNR detections. Second-order bounds have been previ¬ 
ously used for compact binary coalescence waveforms in 
quantifying the accuracy in intrinsic parameters as well 
as the direction of arrival for a network of laser interfer¬ 
ometers [MI]. 

The benefits of using the second order of the expan¬ 
sions is in the fact that they depend up to the fourth 
derivative of the likelihood function and, therefore, are 
sensitive to asymmetries and side lobes of the estimator 
probability distribution (similar to the change in the ac¬ 
curacy of a Taylor expansion when extended to higher 
orders). Also, in the past [5HTT]. the comparison of the 
second order with the first order provided an analytical 
understanding of the reasons the GRLB could not be met 
(for example, in Ref. [10) . a novel relationship between 
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the Kurtosis of the probability distribution of the esti¬ 
mator and the SNR was derived to understand when the 
CRLB could be met). 

Bayesian methods were recently applied to test modi¬ 
fied GR signals through consistency tests [13 [m, and the 
ppE framework [18]. Refs. [muz] developed a framework 
to detect GR violations without modeling the violation, 
this works in the limit of large number of detections. This 
framework was used in GR tests from the GW150914 
transient [6|. Bayesian selection methods were also used 
in Ref. [TH] and Ref. [T5j to constrain the range of ppE 
parameter values, provided that priors are adopted. 

When Bayesian uncertainties are smaller than the fre- 
quentist bounds, it means that the parameter estimation 
errors depend critically on the priors. This issue can be 
an artifact if the prior is not based on previous detections 
or no robustness studies were performed with respect to 
the choice of the priors (see for example the discussion in 
Ref. [5D| of the effects of priors). In this paper, we show 
that this instance happens for an equal-mass binary black 
hole system in the massive graviton case. This example 
illustrates how the present work provides a unique under¬ 
standing of the parameter estimation errors. Although 
GW150914 had a SNR~ 24, its inspiral stage falls within 
the prescribed study of SNR< 20. Q 

In addition, this work extends the Fisher informa¬ 
tion based results of Ref. [IMD, which perform error 
estimations by modifying PN coefficients. We also ex¬ 
tend Fisher-based assessments of specific alternative the¬ 
ories [2TVI25] . Specifically, this paper considers phase 
modification in the restricted ppE framework [S], con¬ 
sidering the ppE framework as a general enhancement 
to existing TaylorF2 [261 W\ GR templates in a three 
detector LIGO-Virgo network m- Galculations in this 
limit were chosen since deformations to the GW’s phase 
are expected to be more resolvable [ia[iHi and comple¬ 
ments recent Bayesian methods testing deviations from 
GR [m [T7]. Second-order frequentist constraints pro¬ 
duced in this paper are at the same order of magnitude as 
the Bayesian model selection’s errors in Ref. ng, where 
our errors are quantified at the one sigma level. As error 
estimates of ppE parameters grow, second-order errors of 
parameters such as the chirp mass, symmetric mass ra¬ 
tio, and time of coalescence also inflate. The results pre¬ 
sented here, and the rescaled bonds which can be simply 
derived by changing the SNR, will be important bench¬ 
marks for any parameter estimation scheme which will be 
used in existing and future interferometer data, including 
Bayesian parameter estimation algorithms. 

Section jl A| of this paper introduces the signal model 
used. Section |ll| discusses what is the resolvable parame¬ 
ter space and the expansion model, in particular subsec¬ 
tion jll^ discusses alternative theories of gravity covered 
in this paper and the asymptotic expansion of the max- 
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imum likelihood estimator model is discussed in subsec¬ 
tion II B| Finally, section O assesses the results, as ap¬ 
plied to a two-dimensional ppE parameter space (IIIA) 
and a seven-dimensional parameter space of equal mass 


(III BI and unequal mass (III GI systems with physical 
parameters included. Results are applied to existing al¬ 
ternative theories of gravity in |III D[ including massive 
graviton, Brans-Dicke, and quadratic modified gravity 
(encompassing Einstein-Dilation-Gauss-Bonnet gravity). 
A summary and discussion is given in section [Iv| 


A. Signal Model 


The waveforms are assumed to be produced by a non¬ 
spinning binary system with all orbital eccentricity in¬ 
formation lost when entering the frequency bandwidth 
of Adv. LIGO and Adv. Virgo. Fourier transform of the 
signal, through stationary phase [13 ESI, becomes, 

4a(/) = , /</„e.« (1.1) 


for the inspiral stage of the compact binaries. For 
the phase '0gr(/) and amplitude Aqj^(/) the standard 
TaylorF2 model [13 [17] is used. 

The signal of a collection of alternative theories of grav¬ 
ity is modelled as (1.1) modulated in the phase and am¬ 
plitude as: 


^gr(/) ^gr(/) (1 +'^^(/)) I (1-2) 

V’gr (/)^ (/) + 


where SA{f) and Si/jif) are a general series of scaling 
parameters ai^Pi € St and in some instances arguments 
call for integer exponentials of i/ry 1/5 [aiEi], where ly = 
(7rM/)i/5 for total mass M and r] = mim 2 /■ Here the 
analysis is done at leading order in the ppE parameters, 

^^ppe(/) = a(j^??i/5)“, (1.3) 


At each interferometer the signal is assumed to be 
recorded with additive noise as in Ref. m- Frequency 
dependent noise for Adv. LIGO are interpolated from the 
official power spectral density [53] of high-power, zero¬ 
detuning. Adv. Virgo is assumed to have the sensitivity 
given in Ref. [34] . For error analysis, and upcoming in¬ 
tegrations, the lower cutoff frequency is set to and 
the upper cutoff is set to the upper limit for reliability in 
the inspiral of the waveform template, i.e., the innermost 
stable circular orbit (ISCO) frequency. 


= 20 Hz , = /:sco « (63/2 ^M)-i. 

The convention used in (I.I) is presented in Appendix [A| 


II. PARAMETER SPACE AND EXPANSION 

For non-spinning systems thirteen parameters are nec¬ 
essary in the description of the inspiral of two coalescing 
















3 


binaries: two mass terms, four angles (two source loca¬ 
tion and two waveform angles), two coalescence parame¬ 
ters, distance to the source, and four ppE parameters in 
the leading order approximation. Singular Fisher matri¬ 
ces might appear [iniiis], indicating that the resolvable 
parameter space is smaller (where the Fisher matrix ap¬ 
proach can still be used). 

The distance is excluded from the error estimates 
because the amplitude has a dependency on both mass 
and distance parameters, and the independent treatment 
of both is unresolvable as already indicated in Ref. m- 
The coalescence phase is also not included because es¬ 
timations of (j)c is relevant only when a full waveform 
(inspiral, merger, and ringdown) is implemented. The 
polarization ij) is excluded because results tend to be in¬ 
dependent of it m- 

Derivatives of the htting factor {FF) [55], 

. [ , \ 

< V\/<«(A)i«(A))\/(«(ai»2(c))/ 

with respect to the binary’s inclination e evaluated at, 
or in a neighborhood of, e = 0 are roughly zero leading 
to impossibility to estimate e and singular Fisher matri¬ 
ces. Here the (-j-) represent noise weighted inner prod¬ 
ucts HUES] and si .2 are GW signals controlled by general 
parameter space vectors A and C- Keeping other parame¬ 
ters fixed and varying only e produces change in the SNR 
equivalent to the rescaling of the distance, which affects 
GW plus-cross polarizations similarly. Top panel in fig¬ 
ure shows the sky-averaged SNR plotted as a function 
of inclination e (only the GR polarizations are consid¬ 
ered). Also, sky patterns of the errors remain consistent 
when varying e. Therefore, since e is degenerate with 
it is also excluded from our resolvable parameter space, 
which becomes = {r 7 ,logAd,tc,lat,long}. 

Throughout this paper amplitude modulations are to 
be held fixed to that of GR: a = 0, because the same 
effect could be produced by changing physical parame¬ 
ters like distance or mass. Such an approach supposes 
that GR-violating amplitudes in the waveform are sup¬ 
pressed or modifications manifest only in waveform prop¬ 
agation!^ Also, recent work suggests that GR modifi¬ 
cations produced during the generation of a waveform 
can be disentangled from that produced during propaga¬ 
tion |31j . thus, in the event that phase deformation dom¬ 
inates GR-violating effects, amplitude modihcations can 
be disregarded. Galculations in this restricted framework 
are performed with modifications at various PN-orders in 


^ Modifications to just propagation could surface through alter¬ 
ations in the dispersion of the GW, with alterations stem¬ 
ming from waveform generation excluded 12211251 . Past studies 
also indicate modulations are most sensitive to phase modula¬ 
tions [T11[28]. 
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FIG. 1. Top: Sky-averaged SNR plotted with e varied for 
system parameters: mi = m 2 = 10Mq,G = (pa = 0, = 

—0.2, Dl = 1100 Mpc, and b = —3 in the three detector 
network. Bottom: Fitting factors (2.11 for a range of /3 with b 
fixed to produce PN-order 0.0, 1.0, and 1.5 modifications for 
a system of: mi = m 2 = IOMq and ta = (j>a = 0. Adv. LIGO 
noise is assumed. Since the range of /3-values scale differently 
at each PN-order, each /3-interval is scaled (as labeled in the 
legend). For example, in the PN-order 0.0 modification the /3 
values in the domain are each scaled by 10“^. 


the phase, where in the strong-field regime discrete val¬ 
ues of b controls what PN-order correction is constituted 
for free parameter j3 (GR result: /3 = 0). 

A qualitative way to study the influence of ppF pa¬ 
rameters (/3, h) on a GR signal can be obtained through 
the correlation of the signals by means of the fitting fac¬ 
tor (13. Fach integration is done from 20 Hz to /isco 
with the noise curve of Adv. LIGO [35] “high-power- 
zero-detuning.” Our exact waveform si is represented 
by a TaylorF2 w avefo rm, w here as, a modified-TaylorF2, 
formed through (1.2) and (1.3), acts as S 2 . So A is the 
GR-limit parameter space vector and C is that of the 
ppE parameter space. The inner products are maximized 


over evenly spaced parameters C to provide a FF-value, 
where FF = 1 represents an exact match between sig¬ 
nals. Both TaylorF2 models are ke pt to PN-order 3.5 in 
the phase. In the denominator of ( |2.1| ), amplitude pa¬ 
rameters normalize to leave jSh in each integrand. 
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The numerator retains integrand (/ 
where, 


A'0(/; X C) = i’{f\ A) -'(/'(/; C) - 

and, in fixing b and varying /?, the parameters needing to 
be maximized over are C = {tc,4>c,Vi ^tat}- Parameters 
are evenly spaced, in a 30 x 30 x 30 x 30 grid, within inter¬ 
vals: 0.05 < r] < 0.25, 0.5Mtot < < 1.5Mtot ,—tt < 

(pc < TT, and —1.3 x 10“^ <tc< 1.3 x 10“^. 

Figure [2 displays the results for an equal-mass system 
of mi = 7712 = IOMq and ta = pa = ^ for PN-order 0.0, 
1.0, and 1.5 modifications in the waveform. Parameters 
C are maximized over for a variety of /3-values. Note that 
at lower PN-orders the interval of /3 is scaled differently 
than the —5 < /3 < 5 depicted, an interval valid for 
PN-order 1.5 modifications. The general trend is that 
the fitting factor is less affected by /3 for larger PN-order 
with a skew in the F'F’-distribution towards the positive 
domain of /3-values. 


A. Restricted ppE template and existing dephasing 
alternatives 

As stated, variations of /3 are restricted to fixed PN- 
order corrections in the phase. For the two-dimensional 
study b is fixed to induce modifications at (separately) 
PN-orders 0.0, 0.5, 1.0, 1.5, 2.0, and 3.0 which acts 
as a demonstration to the error estimation procedure. 
Higher-dimensional studies specifically target a PN-order 
1.0 modification and a weak-field b = —7 modification to 
address dispersion modification and dipole gravitational 
radiation. From this reason /3 is varied with error estima¬ 
tions performed at each /3-value. In Ref. m an analysis 
of binary pulsar PSR J0737-3039 [3H] placed bounds on 
ppE parameters (for this binary 4?/ « 1 as determined 
from radio pulsar measurements [38]). At PN-order 2.5 
(6 = 0 ) degeneracies occur with other fiducial parame¬ 
ters, thus is not considered in the analysis. In some the¬ 
ories constraints for b = —7 cannot be implemented from 
pulsar measurements, due to /3’s dependence on mass dif¬ 
ferences of the system and other theoretical parameters 
which will be discussed shortly. With the exception of 
b = —7, parameters that probe weak-field (6 < —5) are 
not considered since they are better constrained via bi¬ 
nary pulsar measurements [18] . 

At 6 = —7, the even-parity sector of quadratic 
modified gravity (QMG), an example being Einstein- 
Dilation-Gauss-Bonnet (EDGE) gravity, can be explored. 
For even-parity QMG, the violating term for a BBH 
system depends on the mass differences of the BHs: 
/3 oc ~ 47 /), unresolvable for equal-mass sys¬ 

tems 132] . For BHNS systems, the violating coeffi¬ 
cients depend on the ratio of the two bodies: f3 oc 
C 3 ^~*^^(™ns/wbh)^ due to the ‘scalar charge’ vanishing 
in NSs [SniSn]- With this same b = —7 correction, ex¬ 
amples of dipole gravitational radiation, like Brans-Dicke 


(BD), can also be assessed. Here BD-like modifications 
further depend on the difference of parameters which 
measure the body’s inertial mass variations with respect 
to the local background value of the effective gravita¬ 
tional constant. These so-called ‘sensitivity parameters’ 
Sbh.ns are generally set to 0.5 for black holes, so their dif¬ 
ference vanish for a BBH system. Only a BHNS system 
would allow constraints of BD-like modifications since 


0.2 < Sns < 0.3 

For corrections at b 7 ^ —7, most existing modifying 
coefficients depend on parameters that either vanish in 
the non-spinning model (1.1) or contribute beyond PN- 
order 3.5. This is the case in specific models of QMG, 
e.g., the odd-parity sector and dynamical Ghern-Simons 
(GS) gravity |33|. As an example, in the circular inspiral 
of two comparable mass BHs the GR-deviating term of 
dynamical CS has dependencies on the BH spins 51,2 
and their relations to their orbital angular momentum 
L: 6C = 6C{mi^2, Si,2, L) |3S|. When the binary system 
is non-spinning, modifications are beyond PN-order 3.5. 

Beyond modifications during waveform generation, two 
propagating effects are massive graviton (MG) and sim¬ 
plified versions of Lorentz-violating (LV) theories [^I25j . 
Parameters to constrain are the graviton Compton wave¬ 
length Xg and Alv = 27rA^/fo“^). Here A is a phenomeno¬ 
logical parameter modifying the gravitational waveform’s 
dispersion relation. The 7 -dependent distance measure 
D-y (see Ref. [3S] for exact formula) further depends 
on known astrophysical parameters (Hubble parameter, 
matter density parameter, etc.) which are assumed to be 
exact knowns in the analysis |46j . Parameter 7 governs 
the order of correction and 7 = 0 (PN-order 1.0) is what 
we’re limited to since this is the only value contained in 
the ppE framework for the PN-order 3.5 TaylorF2 model. 
Such MG-LV interpretations are generic models modify¬ 
ing the dispersion of a GW with more specific generation 
mechanism still yet to be explored. Ref. [H] notes some 
limitations in prescribing MG effects as modifications of 
the dispersion of the waveform. In LV-type modification 
further work in existing, model-independent approaches, 
e.g., the Standard Model Extension gZl EE), could be 
interesting (see for example Ref. [in])- 

Constraints have been imposed on the wavelength of 
the graviton. The detection of GW150914 and binary- 
pulsar constraint serve as dynamical bounds while solar- 
system constraints, serving as static bounds, provide the 
most reliable estimates PfSDj. So, parameters are repre¬ 
sented by. 


( 10^3[kjn], dynamic (GW), 

Alv = 27rA“^/^, Ag > < 1.6 x 10 ^°[km], dynamic (pulsars), 
1^2.8 X 10^^[km], static. 

For EDGB gravity, the constraint parameter is jaEDCBl- 
Here (3 = , with /3 bbh oc 

^3^-i8/5(i _ 4 ^) and P 

BHNS oc pzr] ®'^®(tons/wbh)^- In 

Brans-Dicke theory /3 oc (sbh.ns ~ 'Sbh,ns)^Wbd ■ From 
measurements of the Cassini spacecraft Emsn bounds 
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on EDGE and Brans-Dicke parameters are, 
laEDGBr^" < 8.9 X 10® km, 

Wbd > 4 X 10^. 

With other suggested constraints [101 [S3] giving, 

laEDGel^^^ < 1-9 km. 
laEDGE^^^ < 9.8 km, 

GW150914 results have allowed studies to infer the theo¬ 
retical signihcance of the testing GR study [S] in various 
specific models, see for example Refs. [SdlISS] . 


B. Asymptotic Expansions 


Similar to Ref. HD, we reasonably assume only Gaus¬ 
sian noise at time of the signal and that the noise is un¬ 
correlated at different interferometers. Here we use the 
analytic asymptotic expansion of the variance and bias 
developed in Refs. [SHII], 

4,=a%ll]+442]+---, ( 2 . 2 ) 

b^i = b^i [1] -I- [2] -I- • • • , (2.3) 

with being the diagonal elements of the covariance 
matrix, where 

for network SNR p. This inverse proportionality contin¬ 
ues at higher orders in similar fashion. Here the network 
SNR is the sum over the square of the optimal SNR 
of the signal at the /-th detector, 

/ 


Notice that p increases for a fixed source by increasing 
the number of detectors. The first-order term of the ex¬ 
pansion of the variance, the diagonal components of the 
inverse Fisher matrix, dominates the bound on the er¬ 
ror in the limit of large SNR, while higher order terms 
become more important for medium to low SNR. 

What is usually regarded as the error in a lab measure¬ 
ment is the square root of the mean-squared error (MSE), 
where the MSE is the sum of the variance (2.2) and 
square of the bias (2.3): MSE^i = a'^i + 6^;. Since this 
analysis computes errors at second-order of 1/p, the ex¬ 
pression above only requires first-order of the bias which 
is negligible as already discussed in Ref. m- We esti¬ 
mate uncertainties of the two-dimensional ppE parame¬ 
ter space for different /3 at a hxed exponential b. In 
addition, the inclusion of 0*^^, to a signal’s extrinsic and 
intrinsic parameter space assessed. 

Finally, error bounds are indicated with. 


Ai9i[l] = Adi [2] = \J [2] 


(2.5) 


For example first-order errors of the symmetric mass ra¬ 
tio p are marked by Ap[l], second-orders are marked by 
Ap[2], and total error with the inclusion of second-order 
contributions as Ap[l -|- 2]. 


III. RESULTS 

In this section we explore the error bounds both as a 
function of the SNR and sky location of the source. The 
asymptotic expansion approach is first applied to a two- 
dimensional ppE parameter space (when the physical pa¬ 
rameters are known) of equal-mass systems. Only phase 
corrections are assumed through unknown ppE param¬ 
eters (/3, 6), while b probes modifications at PN-orders 
0.0-3.0 of the TaylorF2 model (of a PN-order 3.5 phase). 
Based on Ref. IHII] this approach is expected to give 
overly optimistic errors. The Fisher information error 
estimates presented here for the ppE parameters are at 
least an order of magnitude smaller than results with 
Bayesian model selection [15] . 

To identify SNR dependencies and regions of lowest er¬ 
ror estimates the sky dependencies of errors are observed 
through a 289-point sky grid. A point (lati,longj) in 
latitude-longitude coordinates (of the Earth frame) on 
the sky grid follows from the procedure of Ref. E] (de¬ 
tector coordinates also follow Ref. which are fixed 
in the Earth Frame as given in Ref. [561157]). 

As discussed in Section |Hj e = 7r/6 is a fixed value and 
excluded in error analysis. Parameter ip is also fixed and 
arbitrary values can be chosen for fiducial parameters 
(pc and tc- The sky-averaged SNR is restricted to an 
inspiral phase p < 20 to focus on the more likely advanced 
interferometer scenarios. For each system considered, the 
distance of the resolved signal in the network is varied 
to keep a fixed SNR. For a three-detector network (7 = 
77, L, V) the following is chosen for the equal-mass binary 
systems: 

• BBH 1:1- (mi,m 2 ) = (10,10)Mq, Dl = llOOMpc, 

• BNS- (mi,m 2 ) = (1.4,1.4)Mq, Dl = 200Mpc. 

Here the constructed binary black hole (BBH) and binary 
neutron star (BNS) system leaves the network with an 
averaged SNR of p = 14.6 and p = 17.0, respectively. 
For unequal mass systems we choose a BBH system with 
a 1:2 mass ratio and a black hole neutron star (BHNS) 
binary with the following: 

• BBH 1:2- (mi,m 2 ) = (5,1O)M0, Dl = 850Mpc. 

• BHNS- (mi,m 2 ) = (1.4,1O)M0, Dl = 450Mpc. 

which respectively give SNRs of p = 14.9 and p = 15.8. 
For direction reconstruction and related extrinsic param¬ 
eters the network geometry is important; however, for 
intrinsic parameters (as with the ppE parameters) SNR 
gains and losses have a larger impact m- 
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FIG. 2. Sky-averaged errors as a function of P for a two-dimensional ppE parameter space for the BBH 1:1 system of averaged 
network SNR p = 14.6. SNR results of p = 29.3 are also showed by setting the distance to Dl = 550 Mpc. As noted in Ref. m 
error estimates are rescaled as ct[1](p*/p) and cr[2](p*/p) , where p* is the SNR that error estimates are originally calculated 
from. In the top panel the far left column represents each system for a PN-order 0.0 modification {b = —5), the center column 
is a PN-order 0.5 modihcation {b = —4), and far right column is for PN-order 1.0 modifications {b = —3). Similarly, the bottom 
panels are resulting modifications at PN-order 1.5 (6 = —2), 2.0 (6 = —1), and 3.0 (6 = -t-1). P is more tightly constrained at 
lower PN-orders and the inclusion of second-order errors for {P, b) drastically diverge from Fisher estimates as P ^ 0. 
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FIG. 3. Sky-averaged errors, similar to figure]^ for a BNS system of averaged SNR p — 17.0. 


Error Bounds (System) 

PN-order 0.0 

PN-order 0.5 

PN-order 1.0 

PN-order 1.5 

PN-order 2.0 

PN-order 3.0 

A,d[l] (BBH 1:1) 

2.70 X 10"^ 

1.36 X 10“^ 

6.59 X 10“® 

3.07 X 10"^ 

1.39 X 10"^ 

2.66 

A/3[l] (BNS) 

1.29 X 10“® 

1.24 X 10“^ 

1.14 X 10“® 

9.78 X 10"^ 

7.93 X 10"^ 

4.49 


TABLE I. Constant slopes of first-order error bound estimates of the BBH 1:1 (for SNR p = 14.6) and BHNS systems for all 
P values. Here percent errors [%] follow a 1/P relationship for A/3[l] represented above for respective PN-orders. 


In the seven-dimensional study, /3 is varied along b = 
—3, —7 for a BBH 1:1, 1:2, and BHNS systems. The rea¬ 
son for b = —3 is that it simulates modihcations to the 
dispersion of a GW (e.g., massive gravitons or Lorentz vi¬ 
olations [221I25]). Also, b = —7 simulates weak-field mod¬ 
ihcations for dipole gravitational radiation (e.g., Brans- 
Dicke [HUH]) and the non-spinning, even-parity sector 
of quadratic modihed gravity (e.g., Einstein-Dilation- 
Gauss-Bonnett, or EDGB, gravity m)- Distinguisha- 
bility from GR is denoted as the condition that errors 
are smaller than the separation between parameters of 
the GR-limit and that of some alternative theory. 


A. Two-dimenstional study: equal mass 

In this subsection uncertainties for a two-dimensional 
parameter space are computed for both the BBH 1:1 and 
BNS systems, marked by Parameter b is chosen 

at a hxed PN-order correction with PN-order 0.0, 0.5, 
1.0, 1.5, 2.0, and 3.0 (i.e., b = -5,-4,-3,-2,-1,-hi) 
while P is varied at each PN order. Here P probes values 
small enough to induce a sky-averaged error larger than 
100% in b and large enough for < 10% sky-averaged error 


in p. Errors for the BBH 1:1 system are depicted in 
hgure each labeled column representing a particular 
PN-order modihcation. Furthermore, to demonstrate the 
SNR dependence the BBH 1:1 system contains values 
for the scenario in which the SNR is doubled, for this 
the distance is decreased to = 550 Mpc. Figure 
illustrates similar results for the BNS system. 

The constant slopes of errors at hrst-order are cata¬ 
logued in Table |T] for each PN-order. The computed first- 
order errors are consistent with statements of Ref. [TB] 
which demonstrate that different PN-order corrections 
lead to different feasible constraints on /3-values. BNS 
systems offer tighter constraints on P at each chosen b. It 
is interesting to observe that scaling parameters control¬ 
ling propagating modifications, e.g. the graviton wave¬ 
length /3 mg oc are not more tightly constrained with 
BNS systems at shorter distances than BBH systems at 
larger distances. Rather, parameters like /3 mg, also de¬ 
pend on a distance measure and masses of the compact 
objects that adversely affect constraints at shorter dis¬ 
tances and smaller masses. 

The smaller /3, the more second-order effects in the er¬ 
rors contribute. Second-order effects on the errors of b 
are less significant, and only errors > 100% on /3 force 
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FIG. 4. Sky-averaged uncertainties for the equal-mass BBH 1:1 system for a PN-order 1.0 modification of the seven dimensional 
parameter space (ppE parameters {P,b} and physical parameters { 77 , log Ad, G, lat, long}). In the left column the top panel 
displays A/3 percent errors as a function of j3 (the sign of /3 provides different error estimates) and below that are A 6 errors as 
a function of /3 (the sign of /3 does not play a role in these error estimates). In the middle and to the right are the physical 
parameters’ errors, where the constraint of /3 primarily affects the second-order contributions. Enlarging the parameter space 
increases error estimates from those computed in figure at PN-order 1.0, thus weakening constraints on /3. For negative (J, 
the full-dimensional study states A/3[l] = 100% at /3 = —0.16 and A/3[l -I- 2] = 100% at /3 = —0.32. 


sizeable second-order contributions in b. If b is near dis¬ 
tinguishable, A5[l-|-2] < 100%, A/3[l-|-2] are much larger 
than A/3[l]. Only when A&[1 + 2] < 10% do A/3[l] and 
A/3[l-|-2] converge to similar estimates. Simulations pro¬ 
ducing the results of figures and used both ±(3 values 
and the skewed representation of hgure[^is not apparent. 
Note that the range of /3 values, in which error bounds 
are < 100% (figures and [^, are orders of magnitude 
smaller than the ,5-value ranges considered in previous 
studies based on Bayesian methods m- 


B. Full parameter space: equal mass 

The most realistic results come from the study of the 
largest resolvable parameter space. In this subsection, 
first- and second-order uncertainties Ad* of a full 7- 
dimensional parameter space are calculated for the equal- 
mass BBH 1:1 system, where d = {dppE, dphya}- Here b is 
fixed to induce a PN-order 1.0 modification (6 = —3). 
Such corrections simulate effects produced by modify¬ 
ing the GW dispersion relation 0115]. Unlike the two- 
dimensional cases, the errors (first- and second-order) are 
effected by the sign of /3, where sky-averaged errors for 
the ppE parameter pair (/3, b) are displayed in the left col¬ 
umn of hgure Errors of physical parameters affected 


by varying /3 are depicted in the middle and right col¬ 
umn of hgure The skewed behavior of ±/3 results are 
representative of htting factor results of hgure 

For P the hrst-order errors are not at a constant slope. 
A/3[l] approximately follows linear relationship: A/3[l] « 
0.046|/3| -b 0.15, for negative /3. Here a 100% threshold 
error occurs at /? = —0.16, for A/3[l], and /3 = —0.32, for 
A/3[l -1-2]. In this more realistic scenario, it can be seen 
that for extremely small /3 values b falls within its own 
uncertainty. Yet, analogous to the two parameter space, 
a 100% error in A6[l -I- 2] requires large errors in A/3[l -b 
2]. Furthermore, error estimates are at least an order 
of magnitude larger. Another aspect of considering a 
full-dimensional parameter space are the additional error 
trends imparted on physical parameters (masses, arrival 
time, etc) when /3 is varied, see the middle and right 
column of hgure 

The sky distributions of the errors and the SNR are 
shown in hgure 0 Table |ll] catalogs this for —/3 = 
0.25, 0.35, 0.55. This SNR dependence is similar to intrin¬ 
sic parameters for GWs. The /3 values, being a PN-order 
1.0 correction characterizing massive graviton dispersion 
tests, are chosen for the following reasons: 

1. At /3 = —0.25, hgure identihes the conditions: 

A6[2]/A6[l] « 1 with A/3[l] < 100% < A/3[l -b 2]. 

Sky averages are performed before computing the 
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FIG. 5. Sky distribution of error estimates. Color bars represent range of ppE quantities labeled (a), (b)..., (f) in Table [n| 
This demonstrates the correlation of the SNR and ppE error estimation over the sky. See text for discussion. 


ratios. In SNR > 15, we have A6[2]/A6[l] < 1, as 
seen in (a), (b) diplays A/3[l + 2], which ranges 
from 66.4% to 468.7%. A/3[2] dominates the error 
budget. 


ppE /3-value 

Error Estimations 

= 20.8 

pmin — 7.0 

-0.25 




(a) 

A6[2]/A6[l] 

0.55 

1.67 


A6[ll 

12.1 [%] 

36.2 [%] 


A6[l + 2] 

13.8 [%] 

70.5 [%] 


Ab[2l/A/3[1] 

1.19 

3.57 


A/3[l] 

42.7 [%] 

126.4 [%] 

(b) 

A^[l -b 2] 

66.4 [%] 

468.7 [%] 

-0.35 





A6[2]/A6[l] 

0.43 

1.28 


A6[ll 

8.7 [%] 

25.8 [%] 


A6[l + 2] 

9.4 [%] 

42.0 [%] 


A/3[2]/Ad[l] 

0.91 

2.72 

(c) 

A/3[l] 

31.4 [%] 

92.9 [%] 

(d) 

A^[l -b 2] 

42.4 [%] 

269.1 [%] 

-0.55 





A6[2]/A6[l] 

0.32 

0.99 


A6[ll 

5.5 [%] 

16.4 [%] 


A6[l + 2] 

5.8 [%] 

23.2 [%] 

(e) 

A/3[2]/Ad[l] 

0.65 

1.96 


A/3[l] 

21.1 [%] 

62.4 [%] 

(f) 

Ad[l + 2] 

25.2 [%] 

137.3 [%] 


TABLE II. Maxima and minima of estimates depicted in the 
sky-map plot (figure]^ for respective /I-values of figure]^ Er¬ 
rors are the smallest for p„ax = 20.8 and largest for = 7.0. 
Terms labeled with (a), (b)..., (f) correspond to respective 
color bars in figure Values are chosen because they offer 
the most insight. 


2. For (3 = —0.35, sky-averaged A/3[l] < A/3[l -|- 2] « 
100%. Although A6[2]/A6[l] > 1, in limited por¬ 
tions of the sky, the ratio never exceeds 1.3 with a 
maximum of A6[l -|- 2] = 42.0%. There is a strong 
increase in A/3[l-|-2] from A/3[l] in low SNRs. The 
majority of the sky is dominated by second-order 
terms, with A/3[2]/A/3[l] ranging from 0.91 to 2.72. 

3. /? = —0.55 is where we calculate sky-averaged ra¬ 
tio A/3[2]/A/3[1] « 1 with A/3[l] < A/3[l -I- 2] < 
100%. Here larger portion of the sky has ratio 
A/3[2]/A/3[l] < 1 as shown in (e). A majority (but 
not all) of the sky-map has total error falling be¬ 
low 100% after the inclusion of second-orders with 
sky-averaged error at A/3[l -h 2] « 47%. 

From the known dependence on p, quantities displayed in 
figure and Table [H] can be easily re-derived for higher 
or lower SNRs. 


C. Full parameter space: unequal mass 

Here first- and second-order uncertainties Ai3* of a full 
seven-dimensional parameter space are calculated for the 
BBH 1:2 and BHNS system. In this case a weak-field 
b = —7 modification is induced, which in our context 
mimics the non-spinning, even-parity sector of quadratic 
modified gravity (QMG) and can include specifics like 
EDGB gravity. Inclusion of QMG modifications is due to 
j5 being resolvable by a non-zero mass differences at this 
PN-order. These modifications manifest through modifi¬ 
cation of the energy flux as /3 cx C3(l ~ 4?/) and the 
BHNS binary can also test examples of dipole gravita¬ 
tional radiation, like Brans-Dicke (BD). 

Error bounds are presented in figure The overall 
trend of this system’s estimates are similar to the results 
of the equal-mass BBH 1:1 of the previous subsection, 
with a few exceptions. The first being that the separation 
between errors A/3[l], A6[l] and A/3[l -f 2], A6[l -|- 2] are 
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FIG. 6. Sky-averaged error estimates for the BBH 1:2 and BHNS system. Left column represent calculations of the ppE 
parameter errors (A/3, Ab) for negative /3-values, center column are the mass errors (A77, AM), and far right are arrival time Ata 
and latitude-longitude (Alat, Along) error estimates. Here latitude-longitude error estimates are not affected by /3 variation, as 
was previously presented in the equal-mass system. This study states that A/3 bbhi;2[1 -I- 2] = 95.2% at /3 bbhi:2 = —1.8 x 10“"^ 
and A/3bhns[ 1 -t- 2 ] = 95.3% at /3bhns — —4.5 x 10 


not as great as with the PN-order 1.0 modification. In 
comparison to the previous subsection, the chirp mass 
errors AM are roughly the same, yet A?; estimates are 
considerably less. Time of arrival errors At a are also 
less and latitude-longitudinal estimates don’t suffer from 
varying /3 at first- and second-order. 

For the BBH 1:2 system sky contours of ppE and 
mass error estimates at, respectively, |/3| = 1.8 x 10“^ 
and 1,51 = 3.0 x 10“"^ are displayed in figures In 
figure the mass error estimates (bottom color bars) 
are plotted since this /3-value produces sky-averaged es¬ 
timate A/3[l -I- 2] < 100%, with second-order effects in 
the mass estimates making notable contributions (see fig¬ 
ure]^. We observe that in such a context second-order 
effects do not dominate the error budget of Arj and AM 
in this sky-grid. In low-SNR regions, Ary[2]/A?7[l] and 
AA4 [2] /AA4 [1] are near unity. In these same low-SNR 
regimes A/3[2]/A/3[l] > 1 and A/3[I -I- 2] > 100%, which 
demonstrates the sky-grid SNR relation to errors accrued 
on physical parameters due to large error estimates of 
ppE parameters. 

Figure]^ also represents a second set of contours gen¬ 
erated for |/3| = 1.8 X 10“'* modifications. Top color 
bars are representative of ppE parameter error estimates 
{Af3, Ab) valid for this choice of /3. Contours are plotted 
at this /3-value since this simulates the condition that 
A/3[l -|- 2] « 100% with A/3[l] < 100%. Again we ob¬ 
serve the volatility in A/3[l + 2] estimates, ranging from 


53% to about 250% while remaining strongly correlated 
to the SNR. One notable feature of this plot is that ra¬ 
tios A6[2]/A&[1] and A/3[2]/A/3[l] are relatively close to 
each other, being approximately equal to each other in 
regions of high-SNR. This is in contrast to the equal- 
mass study of the previous subsection and demonstrates 
the small separation in A/3[l] and A/3[l -I- 2] estimates 
depicted in the left column of figure which allows the 
ratio A6[2]/A5[l] to be comparable to A/3[2]/A/3[l]. Re¬ 
lations between these quantities depicted in figure can 
be compared to the extrema of the equal-mass BBH sys¬ 
tem of PN-order 1.0 modifications catalogued in Table [IT| 
Similar results come from the BHNS system. 

In order to check that the Fisher information matrix 
did not become singular we systematically explored its 
eigenvalues. For example figure shows scenarios in 
which the Fisher matrix becomes singular for the seven 
dimensional study. These values of /3 were avoided in this 
analysis. 

D. Application to explicit alternative theories 

Since the modification considered in subsection IHIBI 
occur at PN-order 1.0 in the phase, an analysis can be 
done from these results for the massive graviton model. 
Progression of sky-averaged errors for A/3[l -I- 2], cal¬ 
culated from negative ,0-values, of figure imposes a 
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FIG. 7. Sky-map error estimates of ppE parameters A/3 and Ab and mass parameters Ar; and AM for the unequal mass 
BBH 1:2 system. The top color bars for ppE parameters are for /3bbhi;2 = —1.8 x 10“^ and the mass parameters below that 
are for /3bbhi:2 = —3.0 x 10“'^ of results in figure]^ The SNR color bar is valid for both error estimates. Sky-average estimates 
provide A/3bbhi: 2[1 -t- 2] = 95.2%, of /3bbhi:2 — —1.8 x 10 "‘j and A/3 BBH1:2 [1 -b 2] = 47.4% at p BBH1:2 — ~3. 0 X lO"'^. 


constraint of |/3mg| < 0.31. Existing constraints are 
|/3mg, =tatic| < 0.37 and |/3mg, gw| < 2.89 x 10“^, based 
on current static and dynamical (from GW150914 event) 
bounds on Aq (see section |11 A) computed from the BBH 
1:1 system at 1100 Mpc. This asymptotic approach 
thus produces an additional 16.2% constraint on exist¬ 
ing static bounds at Icr. When including second-order 
terms in error estimation the constraints on Xg have a 
fractional increase of 30% from the first-order Fisher ma¬ 
trix approach as calculated in this paper. Given these 


Distinguishability constraint (< 

100 % Error) 

Aj.lv > 3.04 X 10^^ km 

(BBH 1:1) 

< 7.17 km 

(BBH 1:2) 

I^edcbI^^^ < 2.69 km 

(BBH 1:2) 

< 9.45 km 

(BHNS) 

Oedgb < 3.55 km 

(BHNS) 

Wbd > 12.7(sns “ 0.5)^ 

(BHNS) 


TABLE 111. Seven-dimensional study of the BBH 1:1, 1:2, 
and BHNS systems with feasible constraints, i.e., computed 
MSE < 100%. The first considers PN-order 1.0 modifications 
and the latter two consider b = —7 modifications. Included 
are the graviton wavelength (or generic Lorentz-violating) 
dispersion modification and non spinning, even-parity sector 
models of QMG (EDGB parameter included). Brans-Dicke 
constraint depends on sensitivity parameter 0.2 < Sns < 0.3. 


results, further constraints on the graviton wavelength 
Xg may be possible, even with second-order error terms 
accounted for in the low-SNR limit of the inspiral stage 
only. From calculated results the sky-averaged feasible 
bounds are displayed in Table m 

Bayesian assessments in the ppE framework of unequal 
mass systems (of 1:2 and 1:3 ratios) with SNR of 20 put 
constraints at Xg > 8.8 x 10^^ km [TB]. Other Bayesian 
studies also conclude that advanced detecters would gen¬ 
erally not favor a MG theory over that of GR when Xg is 
larger than the most stringent static bounds m- From 
the TIGER method implemented in the testing GR anal¬ 
ysis of GW150914, constraints are at Xg > 10^^ km, 
when the full inspiral-merger-ringdown signal is used (to¬ 
tal SNR of p ~ 24) [ 6 ]. In this respect, our errors impart 
a more conservative approach to error estimation that 
still suggest that constraints may still be improved. 


An application of seven-dimensional results presented 
in subsection 1111 Gl for the BBH 1:2 can also be made. 
This b = —7 modification has /3 qmg oc C 3(1 ~ ^ry). In 
this context the constraint parameter is ps = in 

the non-spinning, even-parity sector of QMG, where ^3 = 
167 ra^p ^3 in EDGB gravity [39]. For the BBH 1:2 system 
figure presents A,d[l] = 99.7% at \f3\ = 1.4 x 10“^ 
and A/3[l + 2] = 95.2% at |/3| = 1.8 x lO"'*. These 


computations translate to respective inputs in Table HI 


for ^3 and Oedgb- Strongest suggested constraints have, 
in terms of EDGB parameter, laEDCBT^^ < 1-9 km and 
laEDGE^^^ < 9.8 km [301 In weak-field tests the 
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FIG. 8 . First order errors (left panels) and eigenvalues (center and right panels) of the Fisher matrix when computations are 
extended to the seven dimensional parameter space. 


Cassini spacecraft has provided |aEDGB|i^^ < 8.9 x 10° 
km (i.e., ^ 3 ^"^ < 2.4 x IQi’ km) [HT| . Bayesian results 
estimate ^ 3 '^^ < 11 km (or IctEDCBr^^ 4 km) at an 
SNR of 20 [32] which is quoted in Ref. [5] as ^3 ^ < 20 
km for an SNR of 10. 

Similar application to QMG and EDGE theories can be 
done with results of the BHNS system. These constraints 


are also presented in Table III and are more stringent 
than the BBH 1:2 system. With BHNS systems Brans- 
Dicke can be investigated through oc (si — S 2 )'^ui~-^, 
where constraint parameter is Wbd with Sbh = 0.5 for 
black holes and for neutron stars 0.2 < Sns < 0.3 gD- 
l44] . Figure 1^ results indicate A/3[l] = 95.3% at |/3| = 
4.5x 10“° for the BHNS system. Thus, constraints results 
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in Wbd > 1-14 and Wbd > 0.51 at Sns = 0.2 and Sns = 
0 . 3 , respecitvely. Results of the Cassini spacecraft have 
also established uj^d > 4 x 10 "^ [55]. In Ref. |5T] Fisher 
estimates placed constants of Wbd > 194 for a BHNS 
systems of similar masses. 


IV. CONCLUSION 


In this paper we implement a frequentist asymptotic 
expansion method to estimate error bounds on the set 
of ppE parameters modifying the phase of the inspiral 
part of low-SNR (p ~ 15 — 17) GW transients. Figure 
provides a summary of the main results of this paper. 
The bound on the mean-squared error estimates from 
compact binaries studied is shown. Each mark repre¬ 
sents the boundary of the (/?, 6)-parameter space where 
the minimum mean-squared error estimates are 100%, 
with P values below each 6-value > 100% and therefore 
not resolvable. Previous Bayesian studies correspond to 
the range of exponential ppE parameter: —11 < 6 < 2, 
as compared to the figure summary. The fact that for 
the massive graviton case (6 = —3) our approach here, 
which is a more realistic lower limit of the Cramer-Rao 
Lower Bound for early detections, rules out results that 
were allowed by a Bayesian study [TB], seems to indicate 
the need of a careful evaluation of the role of the priors. 

Results of the higher order asymptotic analysis of the 
frequentist approach to error estimation states that fur¬ 
ther constraints can be imposed on existing non-GR the¬ 
ories with the study of the seven-dimensional parame¬ 
ter space (see Table III|. This approach does not in¬ 
volve the use of priors. Here the graviton wavelength 
can be constrained by an additional 16.2% as compared 
to current static bounds [50]. Yet, these projected con¬ 
straints do not further bound the graviton wavelength 
when compared to Bayesian estimates or values imposed 


by GW150914 [6|. Note that although GW150914 pro¬ 
vides a constraint of \g > 10^^ km, our result holds for 
a lower SNR of the inspiral stage only. Further stud¬ 
ies present the scenario for the weak-field b = —7 mod¬ 
ification, which can include quadratic modified gravity 
(QMG) (specifics being EDGB gravity) and Brans-Dicke 
type modifications (figure]^. For the non-spinning, even- 
party sector of QMG, bounds suggest further constraints 
are possible as compared to current bounds placed by 
Bayesian estimates and Gassini constraints. Further¬ 
more, error estimates for modifications at both PN-order 
1.0 and the 6 = —7 weak-field follow similar sky-map 
contours, which are correlated to the SNR patterns (see 
figures and . 

General results show that for successively higher 
PN-order modifications, set by 6, the separation be¬ 
tween first- and second-order errors increase (see hg- 
ures|^and[^. Such an effect percolates to the seven¬ 
dimensional study. Error bounds also increase as the 
parameter space is enlarged, where the two-dimensional 
studies provide overly optimistic error bounds. As con- 



PSR .J0737-3039 constraint 
2-dimensional BBH 1:1 
2-diinensional BNS 
7-dimensional BBH 1:1 
7-dimensional BBH 1:2 


7-dimenslonal BHNS 
Massive Graviton [Solar System] 
Massive Graviton [Bayesian] 
Massive Graviton [GW150914] 
EDGB [Bayesian] 


FIG. 9. Constraints on ppE parameters (/3, b). Alongside fre¬ 
quentist mean-squared error < 100% estimates are constraints 
imposed by Bayesian estimates [TH| , solar system tests [^ , bi¬ 
nary pulsar measurements |.37l I38 | . and GW150914 event [6]. 
Regions below each mark/line are where violations cannot 
be detected based on each respective study. The GR-limit 
is /3 = 0. Our frequentist two-dimensional study considers 
ppE parameter space (/3, 6), while seven-dimensional studies 
includes physical parameters (masses, etc.). See text for dis¬ 
cussion. 


straints on /3 become tighter in the seven-dimensional 
studies, the effects of second-order estimates also accrue 
on physical parameters, namely 77, A4, ta, and latitude- 
longitude parameters (see hgures and [^. Finally, 
SNR increases translate error estimates as discussed in 
Ref. [TT| (figure [^, so all results can be rescaled as a 
function of the SNR. 

Galculations performed in this paper are for single de¬ 
tection scenarios. With multiple detections the presence 
of weak, but consistent, violations could be combined to a 
make a stronger statement about error estimations. Such 
methods to resolve consistent signals were explored in a 
Bayesian framework in Ref. m and it is left for future 
studies in the frequentist framework. Furthermore, as 
waveform models advance, for both the inspiral and ppE 
framework, the application of our maximum likelihood 
estimator asymptotic expansion could be applied to spin¬ 
ning binaries or to waveforms that include the merger and 
ringdown phases. This will add insight into additional 
modified theories mappable into the ppE framework. 
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Appendix A: Notation and network signal 


Masses of each compact body are labeled as mi_2, the 
total mass being M = mi -I- m 2 with 1 / = and 

r] = as the reduced mass frequency and sym¬ 

metric mass ratio, respectively. The usual chirp mass is 
A4 = Geometrized units (G = c = 1) are also 

employed [58]. Terms labeled with 1 indicate a partic¬ 
ular quantity for that /-th detector, e.g., is a signal 
received at some J-th detector, is a detector-dependent 
SNR, etc. Finally, the detectors considered are those for 
Adv. LIGO and Adv. Virgo, so we have I = H,L,V for 
the respective advanced interferometers in Hanford USA, 
Livingston USA, and Cascina Italy. Quantities summed 
over / indicate the total network contribution of that 
term, e.g., network SNR, network Fisher matrix. Apart 
from units employed notation follows that of Ref. [351 ■ 
To discuss some of the terms appearing in ([TtJ : T/ is a 
time lag parameter accounting for the delay in the wave¬ 
form’s propagation from the /-th detector frame (IDF) 
to some fiducial frame (FF)j^with and $q being coef¬ 
ficients that depend on the inclination angle e of the bi¬ 
nary system and the generalized antenna patterns 
of each detector. These are represented by. 


U = n • (r/ - rpp) 


(Al) 


( 2 \ 

-I-cos^ e)^ +(-^xCose) j ,(A2) 


, 2FL cos e 

<i)g = arctan ■ 


J-L{1 + cos^ e) ’ 


(A3) 


with n the direction of travel of the waveform, vj the 
distance to the /-th detector (i.e., the IDF origin), and 
Ypp the distance to the FF origin. Reasons for construc¬ 
tion of a frame of common origin is due to the feasibility 
and efficiency displayed in calculations of quantities in 
particular frames. Notion of a common origin between 
the frames is valid since approximative measure^ allow 
the origins of the coordinate systems to coincide. With 
respect to Ref. m the frames are established as the al¬ 
ready mentioned IDF and FF, with a third frame called 
the wave-frame (WF)j^In producing calculable quanti¬ 
ties the frames are then fixed to values of that in the 
Earth frame (EF). 

Since the origins of the frames coincide transformation 
between the frames is feasible through simple Eulerian 
angles with the usual ZXZ convention |^. From this. 


a set of Euler angles converts a quantity from 

the FF into the WF and another set {a^, /3^, 7^) converts 
from the FF into the IDF through the usual rotation ma¬ 
trices. Here angle '(/' is the polarization angle. A variety 
of relations can be uncovered after defining a few new an¬ 
gles. Let angle pairs ($, 0) and (long, lat) describe the 
sources location in the sky (the former being in spherical 
coordinates and the latter in longitude-latitude coordi¬ 
nates), let (S, (^) be defined from projections of n onto 
the FF’s axis, define angles ,T^) so that they pre¬ 
scribe the location of the /-th detector with respect to 
the FF, and allow angle to span the region between 
the first detector arm (in the IDF) and the local northern 
direction. These relations are summarized as follows: 

0 = $_|=long-| = S+| (A4) 

6* = 7r— 0= — -|-lat = C 

and 

a^ = n^ + l, /3^ = |-T^ 7" = A^ + |. (A5) 

Formulation of Tp ^ into a symmetric-trace-free base 
has been performed, with respect to the Eulerian angle 
dependence, and what surfaces in the frequency repre¬ 
sented signal are the two generalized antenna patterns: 

K = ^(T2,(a^/3^7")+7^-2s(a^/3^7")) (A6) 

K = ^(T2,(a^/3^7")+7^-2.(a^/3^7")) (A7) 

where Tmn are second-order Gel’fand functions 
being their complex conjugates). Function statements, 
such as and represent their de¬ 

pendencies on Euler angle rotations from FF ^ IDF 
and FF WF, respectively. See Ref. m for exem¬ 
plary calculations. Note that an auxiliary ppE template 
has been developed that considers extra polarizations of 
waveforms produced in non-GR gravity, incorporating 
additional propagating degrees of freedom in the ppE 
framework |3T| . Although it is of interest to measure ex¬ 
tra polarizations expected in a variety of alternative the¬ 
ories of gravity, these extra modes lead to more complex 
models. For initial analysis of modified gravity through 
the asymptotic maximum likelihood estimator approach 
a ppE template, with only the standard two propagat¬ 
ing modes, is considered both sufficient and satisfactory 
for now. Ref. [60] investigated methods to test non-GR 
polarizations via continuous waveforms from asymmetric 
pulsars. 


® FF is the frame in which the origins are referenced to coincide. 
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